Is the spectrum of highly excited mesons purely coulombian? 
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We show that a static central potential may provide a precise description of highly excited light 
unflavoured mesons. Due to string breaking this potential becomes of chromoelectric type at suffi- 
ciently large quark-antiquark distances giving rise to a coulombian spectrum. The same conclusion 
can be inferred for any other meson sector through a straightforward extension of our analysis. 
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In the last years the interest in the highly excited 
light-quark meson spectrum has been renewed as a con- 
sequence of the observation of more than thirty new me- 
son resonances with masses between 1.9 and 2.4 GeV 
in the exhaustive analyses of Crystal Barrel and PS172 
data [ij . Many of the reported new resonances are listed 
by the Particle Data Group (PDG) in the section 
"Other Light Unflavoured Mesons" awaiting for confir- 
mation from a separate experiment. Regarding the non- 
new resonances also found in [H they are in perfect cor- 
respondence (quite similar masses) with the ones listed 
in the "Light Unflavoured Mesons" section of the PDG. 
From the theoretical point of view the extensive spectrum 
up to 2.4 GeV has made clear an approximate hydro- 
gen like classification of meson states [3|. In this article 
we show through a dynamical study that the physical 
origin of the hydrogen like degeneracy may have to do 
with string breaking once the quark-antiquark average 
distance in a meson reaches a sufficiently high value. We 
shall rely on a Gonstituent Quark Model (GQM) calcu- 
lation by solving the Schrodinger equation for a static 
potential. Although the use of a static potential for light- 
quark systems is in general debatable we shall justify its 
applicability for highly excited states. Then in the spirit 
of CQM calculations we shall assume that relativistic cor- 
rections to the Schrodinger equation may be to a good 
extent taken into account through the effective character 
of the parameters of the potential. 

From lattice QGD the static potential is the ground 
state energy of a bound state composed of two static 
colour sources and gluons 0]. In the so called quenched 
approximation (only valence quark (q) and antiquark (g)) 
the static qq potential, given by the expectation value 
of the Wilson loop operator, resembles a funnel poten- 
tial containing a confining term depending linearly on r, 
the q-q distance, plus a coulombian term proportional to 
(1/r). When including sea quarks, i.e. in the so called un- 
quenched approximation, the static sources q and q are 
screened by light quark-antiquark pairs that pop out of 
the vacuum. Then transitions between the static sources 
state (string) and mesons coming from the recombination 
of the static sources and the members of the pairs may 
take place. Two physical effects occur. On the one hand 
the confining term becomes a constant (string breaking) 
at a certain saturation distance, r^. Although this be- 



haviour has not been detected with the Wilson loop tech- 
nique, the finite temperature potential extracted from 
Poliakov line correlators at temperatures close to the de- 
confinement phase transition, exhibits a flattening, once 
sea quarks are included into the action j5| . Such a flatten- 
ing occurs at distances ~ 1.15 fm.The diagonalization 
of a two by two correlation matrix between the string 
state and a two meson state (each meson formed by one 
static source and one member of a pair) should confirm it. 
On the other hand the coulomb strength remains higher 
than in the quenched case. From these indications we 
shall assume a static potential of the form 
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where a stands for the string tension, k for the coulomb 
strength and C for a constant to fix the origin of the 
potential. 

We shall choose this form for the effective qq poten- 
tial in our GQM calculation. Let us recall that in quark 
models the effective qq potential is found by equating the 
scattering amplitude of free quark and antiquark with 
the potential between bound quark and antiquark in- 
side a meson. Thinking of a single exchange diagram 
the static limit corresponds to no energy transfer, i.e. to 

9° = i^q)initial ~ final ~ final ~ initial ~ ^■ 

So the static approximation means q^ ~ 0. Then by sub- 
stituting Eq — rriqW 1 + {/p'i/raVj we can easily establish 



as a criterium for the validity of such approximation the 
requirement 



initial 



f inal 



« 1 



(2) 



By replacing (p',),„,,,,; ^ ^ + {~P q) f.nai ^he cri- 
terium (m can be written as 



T. << 1 



(3) 



For non-relativistic as well as for relativistic systems ([3]) 
is satisfied when 
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« 1 and 
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Since the main contributions to the interaction come from 
distances r « he/ | g | we expect the mesons satisfying ^ 
to have root mean square radius (rms-radius) 
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Notice that for bottomonium (bb) and charmonium (cc) 
the conditions (O and I'Pql /rriq :^J^ I are well satisfied 
according to reference [6| where a good description of 
such mesons is attained by means of a quenched poten- 
tial. The validity of the quenched approximation in this 
case can be related to the fact that for bb and cc the sat- 
uration distance Tg may be significantly larger than 1.15 
fm. The screening effect on bb and cc has been studied in 
the literature From these studies a value for (''s)f,ft 
up to 1.8 fm may be conjectured. This value may be 
reflecting the fact that the bb and cc decays into stable 
hadrons involve more than two mesons. Then it can be 
explained why the predicted bb and cc spectra involv- 
ing states with < >^^^^ (^s)bhcc hardly change from 
the unquenched to the quenched approximation. The ex- 
perimental extension of the spectra to states with larger 
rms-radii is essential to confirm or refute this conjecture. 

Here we centre on light-quark mesons. We shall re- 
strict for simplicity to isospin 1=1 mesons (they contain 
only u (u) and/or d (d)): ttj, bj, pj and aj (J : total 
angular momentum). Thus we avoid all possible compli- 
cations coming from qq annihilation and from ss com- 
ponents. The mass of the constituent quarks rn^^--^ ^f-^y 
named henceforth m„, is a parameter of our model. We 
shall fix its value from the average dynamical quark mass 
generated by Spontaneous Chiral Symmetry Breaking 
(SCSB), m„(|'g*|), in the energy region under consider- 
ation. From instanton model calculations Q confirmed 
by lattice QCD we know the explicit m^d q |) de- 
pendence so that it has its maximum value at | 5 | = 
(to„(0) ~ 0.350 GcV) and decreases when increasing 
\~q\- For |~^| = 0.1 GcV for instance one has r7i„(0.1 
GeV) = 0.332 GeV. Therefore g]) tells us that the static 
approximation might only be applied for \~q\ < 0.1 GeV. 
Then we shall use m„ = 340 MeV as an average mass in 
this interval for which we expect from ([5|) to have mesons 
with < >i/2>> 0.6 fm. It should be added that SCSB 
has another important effect: the appearance of Gold- 
stone bosons. This effect is not explicitly reflected in ([T]). 
We shall comment on this later on. 

To check the applicability of the static approximation 
we proceed to calculate the spectrum of / = 1 mesons 
with an effective potential of the form ([1]). For this pur- 
pose we have to fix the parameters of the potential. For 
rs we take {rs)u = 1-15 fm. For the string tension ct„ we 
shall use the value = 932.7 MeV/fm obtained from 
the Regge trajectory for pj and aj (see for instance [1]) 
in accord with lattice evaluations. Regarding fc„ and Cu 
we shall fix them by fitting the average masses of the 



experimental states with higher orbital angular momen- 
tum L, since we expect these states to have large rms- 
radii due to the centrifugal barrier. Let us realize that 
the calculated masses coming out from the Schrodinger 
equation will only depend on L {L = 0, 1, 2, 3...) and on 
the radial quantum number Ur {rir = 1, 2, 3...). We shall 
denote them as Mi_„^. So we should compare them to 
the average masses of experimental {L,nr) multiplets, 
{{M)L,nr.)Exp- The existence of such multiplets has been 
suggested elsewhere Actually the assumption of a 
long distance interaction depending only on r drives nat- 
urally to a SU {4) Spin- Isospin X 0(3) group of symmetry 
so that the 1=1 light unfiavoured mesons belong to 
15-plets (note that the product of SU{4) quark and an- 
tiquark representations is 4x4= 15 + 1; the singlet 
representation 1 contains only / = mesons that we do 
not consider). Indeed we should better talk about super 
15-plets since for each member of the multiplet there are 
so many experimental states as possible different J val- 
ues (J degeneracy). As we consider only 1=1 mesons 
we define their average mass in the corresponding 15-plet 
(L, Ur) as 

(Ml,„J£,p= [^(2J+1) j ^(2J+l)Xj (6) 
\ J / x,j 

being Xj the experimental masses assigned to the mul- 
tiplet [Xj = AUj, Mp^ or Afb,,, MaJ. 

The maximum value of L for which we have some can- 
didate from the PDG catalog (the section "Light Un- 
fiavoured Mesons") or from reference [H, named hence- 
forth CBC (for Cristal Barrel Collaboration), is i = 5. 
In fact there is only one candidate, the PDG resonance 
a6(2450±130). Although the error bar is big the existence 
of a non considered 7 = PDG resonance f 6(2465 ± 50) 
that can be assigned to the same 15-plet (L = 5, = 1) 
makes us confident about the average PDG mass. For 
(L = 4, rir = 1) we also have only one PDG resonance, 
Pg(2330 ± 35) but a complete set of CBC candidates 
(7r4(2250±15), P3(2260±20), P4(2230±25), P5(2300±45)) 
(an explanation for the PDG - CBC difference in mass 
for pg is given in reference ll,]). Let us also notice the 
existence of / = CBC resonances, W3, 0^4 and cjg at 
about the same mass. 

By choosing /c„ = 2480 MeV.fm and C„ = 1070 MeV 
we reproduce correctly their average masses (for (L = 4, 
rir = 1) a, value in between the PDG and CBC averages 
is chosen) . It is noteworthy that the calculated rms-radii 
for the fitting states (3.7 fm and 2.6 fm) consistently 
satisfy < >^^^» 0.6 fm. Moreover the calculated 
values of \'p'q\/mq (1.05 and 1.3), although indicating 
the relativistic character of quark and antiquark, are not 
much greater than 1, as required from ([4]). 

The results for the spectrum of states with < >^/^> 
2.0 fm, for which we expect the static approximation may 
work (all the states have | p q\ /triq < 1.7), are shown in 
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Table H] (we include for completeness multiplets with rms 
radii below 2.0 fm as (2,2) and (3,1)). Our results are 
compared to CBC and PDG average masses. The states 
entering in the calculation of the averages are specified. 
The multiplets (1,4) and (1,3) lack an ao despite hav- 
ing an available candidate ao(2025). The reason for not 
including this state is, appart from the ambiguity in its 
assignment, the general deficient description of ao states 
provided by quark models pointing out the need to in- 
corporate more than two valence components. It should 
be noted additionally that /93(2260 ± 20) has been as- 
signed to two multiplets (4,1) and (2,3). The reason 
for this double assignment is the assumption that the 
resonances belonging to such multiplets would be, as 
indicated by their partners in the multiplets, almost de- 
generate. All the multiplets considered have at least one 
PDG cataloged state or one resonance rated at least three 
stars in reference [l|. The multiplet (0,4) has not been 
considered since the only trustable assignment to it, the 
CBC 7r(2070), has only a two-star rating (let us point 
out that if we assigned the PDG /9(1900) to this multi- 
plet our result would be in perfect accord with data). As 
can be checked our results seem to agree with data for 
< >^/^> 2.0 fm. To be more precise we can rely on the 
approximate linearity and equidistance of Regge trajec- 
tories satisfied by data up to a mass of 2.3 GeV [H, [3l.Tlol| 
with a standard Regge slope of about 1.1 GeV^. From 
our results for (1,3), (2,3) a Regge slope of 1.4 GeV^, 
far above the standard value, would be obtained for the 
correponding (i, 3) trajectories. We interpret this as an 
indication that the static approximation is doubtful for 
the (1,3) state. Let us also note that a Regge slope of 0.79 
GeV^, quite below the standard value, would be obtained 
from our chosen masses for (4, 1) and (5, 1) in the (L, 1) 
trajectories indicating that the calculated spectrum tends 
to a coulombian one when increasing the energy. This 
is clearly indicated by the almost mass degeneracy for 
states with the same values of {L -\- rir) for (L + rir) > 6 
(< >i/2> 4 fm). Then we can give a closed formula 
for the mass of highly excited mesons 



{Ml 



T-r ) {L+n,.)>6 



kl 



2 {L + flrf 



Cu 



(7) 



where p — m„/2 is the reduced mass of the system. The 
predicted masses for multiplets with the same value of 
{L + rir) are listed in Table HIl up to (L -|- rir) = 9. Let us 
also remark the accidental degeneracy between positive 
(+) and negative (— ) parity states with L{+) — L{~) = 
odd — nr{—) — nr{+). Moreover, as the coulombian en- 
ergy (third term on the right hand side of ([7|)) has as an 
upper bound we predict that the 7 = 1 light unflavoured 
meson spectrum has an upper bound or limiting mass 



(L, rir) 


^ 2-1/2 

fm 


MeV 


iiM)L,«r)cBC 

MeV 


{{M)L,nr)pDG 

MeV 


(5,1) 


3.7 


2450^ 




2450 ± 130 
as (2450) 


(1,4) 


3.4 


2255 


2219 ±43 
bi(2240) 
ai(2270),a2(2175) 




(2,3) 


3.2 


2254 


2248 ± 37 
7r2(2245),/9(2265) 
Pa (2225), Pa (2260) 


Pa (2250) 


(3,2) 


2.9 


2258 


2258 ± 38 
bg (2245), aa (2255) 
as (2275), a4 (2255) 




(4,1) 


2.6 


2283^ 


2262 ± 28 
7r4 (2250), P3 (2260) 
P4(2230),/95(2300) 


2330 ± 35 
P5(2350) 


(1,3) 


2.1 


1919 


1947 ± 47 
bi(1960) 
ai(1930),a2(1950) 




(2,2) 


1.9 


1913 


1980 ± 23 
7r2(2005),p(2000) 
P2 ( 1940), P3( 1982) 


P3(1990) 


(3,1) 


1.6 


1937 


2023 ± 24 
bs (2032), a2 (2030) 
as (2031), a4 (2005) 


2010 ± 12 
a4(2040) 



TABLE L Calculated masses (column III) and rms-radii (col- 
umn II) for (LjUr) multiplets (column I) from the set of pa- 
rameters niu = 340 MeV, cr„ = 932.7 MeV/fm, fc„ = 2480 
MeV.fm and Cu = 1070 MeV. Experimental CBC and PDG 
average masses (columns IV and V) are shown for comparison. 
In both cases the candidates to be members of the multiplets 
are indicated. The superindex f in the (5,1) and (4,1) calcu- 
lated masses indicates the average mass values chosen in the 
corresponding multiplets to fix the parameters. 



(L + rir) 
(L,n,.) 


mL,r,r{MeV) 


6 

(0,6), (1,5), (2,4), (3,3), (4,2) 


2450 


7 

(0,7), (1,6), (2,5), (3,4) 
(4,3), (5,2), (6,1) 


2548 


8 

(0,8), (1,7), (2,6), (3,5) 
(4,4), (5,3), (6,2), (7,1) 


2613 


9 

(0,9), (1,8), (2,7), (3,6), (4,5) 
(5,4), (6,3), (7,2), (8,1) 


2657 



TABLE II: Predicted masses for some {L,nr) multiplets with 
L + rir > 6. Parameters as in Table IT] 
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given by 

Limit 



+ mu + CT„(r,)„ + C„ = 2823 MeV (8) 



This limit is compatible with current data since reported 
resonances with higher mass listed in the PDG section 
"Other Light Unflavoured Mesons" may be assigned to 
mesons containing ss. Although a specific analysis par- 
allel to the one just performed would be required for 
these states we can expect their limiting mass to increase 
with respect to the value in ([8|) by at least an amount 
[{ms - ruu) + (m^ - rriu)] ^ 300 - 500 MeV. Then the 
resulting hmit (< 3300 MeV) would be compatible with 
all existing light unflavoured meson candidates. Beyond 
this limit one meson states can not exist. Instead the 
system fragments into several mesons. 

Furthermore a mass limit for baryons containing only 
quarks u and d may be derived through the simple pre- 
scription {MB)Lirmt ^ 3m„ -|- (3/2) [(T„ (r^ -|- C„] = 
4213 MeV. This value is consistently above the most mas- 
sive reported nucleon and delta resonances in the PDG 
sections iV(~ 3000 Region) and A(- 3000 Region). 

With respect to the fitted values of the parameters 
some comment is in order. As we do not expect signif- 
icant Goldstone boson contributions for we are dealing 
with q-q distances much larger than the parameters 
may be incorporating mainly gluon contributions (apart 
from possible relativistic quark kinetic corrections). So 
the coulomb strength, k, can be tentatively related to an 
effective quark- antiquark-gluon coupling, a^, through the 
colour relation k = (4/3) a^- Then from the fitted value of 
k = 2480 McV.fm we get asi\'q'^\ < 0.01 GcV^) 9.4. It 
is interesting to realize that this value for is precisely 
the same reported in reference jTlj from a solution of the 
truncated Schwinger Dyson equations. However this "co- 
incidence" should be taken with caution since a smaller 
value has been obtained in other calculations 12l|. As a 



matter of fact the only general conclusion we may extract 
from the work of all the groups is the almost constancy 
of as for sufficiently low | ~q^ | . Undoubtedly, the validity 
of ^ for highly excited light-quark mesons has much to 
do with this constancy of ag for |~^^| < 0.01 GeV^. 

To summarize, our study of highly excited 1 = 1 light 
unffavoured mesons shows that their spectrum could be 
purely coulombian. As a consequence it would have a 
limiting mass. In fact the idea that the meson spectrum 
might have an upper bound is not new. A mass limit of 
3.2 GeV for light-quark mesons involving uu, dd as well 
as ss components was suggested some years ago from a 
study of Regge trajectories ^3]- Our dynamical analy- 
sis supports such suggestion confirming string breaking 
as the underlying physical mechanism. However the lack 
of trustable experimental data beyond 2.5 GeV leaves 
opened the door to other interpretations. For instance 
one could think that the observed flattening of the confin- 
ing potential corresponded indeed to a severe softening of 



the confining interaction in the region under study. Then 
the good effective description achieved would be compat- 
ible with a very slight increase of the interaction with the 
distance and consequently with an unbound meson spec- 
trum. Besides the non-relativistic hydrogen like sym- 
metry we have made dynamically evident may be only 
the effective face of a broader relativistic symmetry in 
the energy region under consideration as discussed in 0] . 
In our model the higher the meson mass the lesser the 
\ p q \ /ruq value and the less relativistic the system (when 
approaching the limiting mass IJi^q] /mq 0) implying a 
non-relativistic coulombian symmetry at the long range. 
Consistent relativistic calculations beyond our effective 
CQM treatment could shed more light on this point. 

Let us also emphasize that the analysis we have per- 
formed may be extended to any other meson sector al- 
though the current lack of data makes it not feasible. 
Then a definite answer to the general question about the 
coulombian nature of the highly excited meson spectrum 
has to be postponed until more complete data are avail- 
able. We encourage an experimental effort along this line. 
In the mean time we hope our results may be suggestive 
and motivate further studies in the field. 
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